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Abstract 

As a sequel to our earlier work on wino-dominated Xi an d X 2 (wino models), we focus on 
the pMSSM models where xf, X 2 an d X 3 are either higgsino dominated (higgsino models) or 
admixtures of significant amount of higgsino and wino components (mixed models), with or 
without light sleptons. The LHC constraints in the trilepton channel are significantly weaker 
even in the presence of light sleptons, especially in the higgsino models, compared to those mostly 
studied by the LHC collaborations with wino-dominated xt an d The modes X 3 JX 2 —t X? ^ 
with large branching ratios (BRs) are more common in the higgsino models and may produce 
spectacular signal in the LHC Run-II. In a variety of higgsino and mixed models we have 
delineated the allowed parameter space due to the LHC constraints, the observed Dark Matter 
(DM) relic density of the universe, which gets contributions from many novel DM producing 
mechanisms i.e., the annihilation/coannihilation processes that lead to the correct range of relic 
density, and the precise measurement of the anomalous magnetic moment of the muon. In 
the higgsino models many new DM producing mechanisms, which are not allowed in the wino 
models, open up. We have also explored the prospects of direct and indirect detection of DM 
in the context of the LUX and IceCube experiments respectively. In an extended model having 
only light gluinos in addition to the electroweak sparticles, the gluinos decay into final states 
with multiple taggable b-jets with very large BRs. As a consequence, the existing ATLAS data 
in the 01 + jets (3b) + fir channel provide the best limit on (« 1.3 TeV). Several novel 
signatures of higgsino models for LHC Run-II and ILC have been identified. 
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1 Introduction 

The first phase of the p-p collision at the Large Hadron Collider (LHC) (Run-I) has given lower limits 
on the masses of the super-particles (sparticles) for models involving supersymmetry [1-3], although 
the latter is yet to be discovered. The bounds on the masses of the strongly interacting sparticles 
are already stringent [4-9] 1 . Searches have just been started at the new round of experiments at 
higher energies (13 TeV). However, the possibility that these sparticles are even beyond the reach of 
the ongoing experiments as well is wide open. It is known that this scenario is indeed favoured by 
the SUSY flavour and SUSY CP problems [2]. On the other hand, we also note that heavy squarks 
belonging to the first two generations do not spoil the naturalness [11,12] of a SUSY model. 

It is remarkable that the observed mass of the Higgs boson at around 125 GeV [13] at CERN is 
well within the MSSM predicted upper limit of Mh (^ 135 GeV), where Mh refers to the mass of the 
CP even neutral lighter Higgs boson h. In this analysis we consider the decoupled Higgs scenario 
of the MSSM, namely Ma » Mz,Mh [2,3], where Ma refers to the mass of the pseudoscalar 
Higgs boson. In the decoupling limit h becomes Standard Model like in its couplings [14]. It is well 
known that the current Higgs data are indeed consistent with the decoupling limit [15]. 

If the heavy squark-gluino scenario along with a decoupled Higgs sector is indeed realized 
in nature, we must accept that the observability of SUSY signals hinges on the properties of 
the sparticles in the electroweak (EW) sector 2 . Although the production cross-sections of these 
sparticles are rather modest, significant bounds on their masses have already been obtained by 
both the ATLAS and CMS collaborations of the LHC [16-21] from the null results of (i) chargino 
(xf) - neutralino (x°) searches 3 via the process pp —l xfx® leading to the trilepton + transverse 
missing energy (fir) signal and (ii) slepton searches via the opposite sign same flavour dilepton 

1 However in compressed SUSY type scenarios limits on sparticle masses are considerably weaker [10]. 

2 The fermionic members of this sector, the charginos and the neutralinos, are referred to as the electroweakinos 

while the scalar members are sleptons of both L and R types and sneutrinos. 

3 Throughout this paper chargino would stand for the lighter chargino (xf) unless otherwise mentioned and the 

four neutralinos x? — xl are arranged in order of ascending masses. 
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+ frp channel. We shall focus on the analyses performed by the ATLAS group. They obtained 
model independent upper bounds on the cross-sections of these processes applicable to any Beyond 
the Standard Model (BSM) scenario corresponding to different signal regions each characterized by 
an appropriate set of selection criteria. The results were interpreted in terms of several simplified 
models. A large number of phenomenological analyses have addressed the electroweakino searches 
and related topics in the context of the LHC [22-24]. 

The discovery potentials of the charginos and the neutralinos depend on their pair production 
cross-sections and decay branching ratios (BRs) into leptonic channels which contribute to the 
trilepton signal 4 . These observables depend - among other things - on their compositions. The 
focus of this paper is on the phenomenology of models where xf, X 2 an d X 3 are either higgsino 
dominated or admixtures of significant amount of higgsino and wino components. In the following, 
we shall refer to : 

i) the former class of models as the xf'-higgsino or simply, the higgsino models , and 

ii) the latter class of models as xf -mixed or simply, the mixed models. 

On the other hand, the LHC collaborations restricted their analyses of chargino-neutralino 
search in the trilepton channel to simplified scenarios where the lighter chargino (x^) an d the second 
lightest neutralino (^ 2 ) are wino dominated and are nearly mass degenerate. All models belonging 
to this class will be referred to as the xf -wino or simply, the wino models. As in the analyses 
of the ATLAS or the CMS collaborations the lightest supersymmetric particle (LSP), i.e., the 
lightest neutralino (xi) is considered to be strongly bino-dominated over the parameter space with 
m~± >> m~ 0 . In the higgsino models this is ensured by the chosen hierarchy /i = M 2/2 >> Mi, 
where Mi, M 2 and /i are the U(l), SU(2) gaugino and higgsino mass parameters respectively. In 
some regions where fj, satisfies n = M 2/2 ^ Mi the LSP gets non-negligible higgsino components. 

In the backdrop of the basic varieties for the compositions of the electroweakinos and the 
correlation of slepton masses with that of the electroweakinos, that we are going to enumerate 
shortly, one of the main goals of this paper is to go beyond the wino models and reinterpret the 
ATLAS data in several higgsino and mixed models. Before we move on to the above models, we 
would like to set the background briefly by digging into the analysis of the xf- wino scenario. 

4 In this work lepton usually implies electrons and muons unless mentioned otherwise. 
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X x - wino: Each simplified wino model considered by ATLAS and CMS analyses [16,19-21] 

belongs to either of the two following broad categories. In one case, sleptons of all three flavours 

are heavier than the winos (the Light Wino and Heavy Slepton model (LWHS)) 5 . In the other 

category, at least one type of slepton of all flavours (i.e., either L or R-type or both) is lighter 

than xt and X. 2 - The BRs of the electroweakinos and consequently, the discovery channels in 

the two scenarios may be significantly different. The second category in turn consists of several 

subcategories depending on the type of the light sleptons, and their masses with respect to m~o 

Xi 

and m~±, some of which were considered by the ATLAS group while the CMS collaboration as well 
as the analysis of Ref. [24] studied more variations. The important features of each subcategory 
will be summarized in the latter sections. In the wino models xt and x 6 production followed by 
their decays into trileptons is the main discovery channel in most cases. 

We wish to stress that in Ref. [24] as well as in this work we do not restrict ourselves to stay 
within the decoupled slepton scenario. The presence of relatively light sleptons has two important 
implications. First, the constraints from the direct slepton searches at the LHC must be included 
in our analysis. Moreover, the possibility that the sleptons play active roles in DM production 6 is 
resurrected, as has already been noted in Ref. [24] and will be further illustrated in this analysis. 

The slepton search results, mainly for the selectron (e) and the smuon (Jl) are fairly insensitive 
to the electroweakino sector. We particularly note that sleptons of the first two generations have 
negligible L-R mixing, which in turn means that there is hardly any dependence of // or tan /3 in 
determining the masses of these sleptons or their couplings with the gauge bosons. Any slepton 
lighter than xt or X 2 decays into its fermionic superpartner and x? with 100% BR and this is 
independent of the composition of xt /Xi- Subject to these general assumptions mass bounds were 
obtained by the LHC collaborations in several simplified models [17,21], In view of the above 
independence we shall directly use the constraints from the slepton searches as derived in Ref. [24] 
for the wino models. 

5 In Reference [24] this model was called the Light Gaugino and Heavy Slepton (LGHS) model but we find the 
terminology LWHS to be better suited for discriminating between the wino and the higgsino models. A similar change 

in nomenclature applies to all the models discussed in Ref. [24]. 

6 The DM producing mechanism which would often be quoted in this work would mean annihilation/coannihilation 

processes that bring the DM relic density within the acceptable range given by the WMAP/PLANCK data. 
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The models studied by the ATLAS collaboration' to interpret their search results are in some 
sense oversimplified. The parameters of such a model provide a minimal set to understand im¬ 
portant aspects of a SUSY signal. However, in order to test SUSY in the light of the LHC as 
well as other constraints from the so called indirect tests, one requires a closely related but a 
complete model like the phenomenological minimal supersymmetric standard model (pMSSM) [25] 
with additional parameters. In Ref. [24] we enlarged each simplified wino model by introducing a 
minimal set of parameters belonging to the EW sector (like g and tan/3) and tried to scrutinize the 
pMSSM thus obtained by considering three major constraints, namely, the LHC mass limits on the 
chargino-neutralino and the slepton sectors, measured dark matter (DM) [26-28] relic density of 
the universe from the WMAP [29]/PLANCK [30] 7 8 and the precisely measured value of (g — 2) M [35]. 
Moreover, we dispensed with the unrealistic assumption m~ = mj which leads to erroneous LHC 
limits especially if m~ o ~ mj . These modifications change the LHC limits quite significantly in 
some models and we computed these changes in Ref. [24] by a PYTHIA (v6.428) [36] based anal¬ 
ysis using ATLAS data. New bounds for several wino-slepton mass hierarchies not considered by 
the ATLAS collaboration were also derived. For each model, compatibility with the three major 
constraints delineates an allowed parameter space (APS). Each APS, in turn, enables us to focus 
on the expected SUSY signals in the future LHC experiments. 

higgsino or mixed : In this paper we extend our earlier analysis to the xf -higgsino 
and xf - mixed models. The mass hierarchies among the sleptons and electroweakinos are, however, 
similar to the ones in [24]. In Table 1 we present the models analyzed in this paper and the choice 
of parameters for each of them. In our analysis, M\ and M 2 are free parameters. 

As mentioned earlier the LSP is also assumed to be bino-dominated with some degree of higgsino 
mixing, depending on the parameter space. Finally, we delineate the APSs in both the models while 
isolating the effects of each major constraint clearly. We emphasize that in the post-LHC era, these 
models, especially the higgsino models , have not received due attention in the literature. Yet, the 
difference in phenomenology of the higgsino models and that of the wino models is indeed worth 
noting. Firstly, the LHC exclusion contours from the trilepton searches shrink significantly in the 
higgsino models and even become irrelevant in some scenarios. As we will see later the physics of 

7 Similar simplified models were also analysed by the CMS collaboration. 

8 A partial list of works on SUSY DM may be seen in Ref. [31-34]. 
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Model 

Acronym 

Parameter Choice 

xf-Higgsino 

Light Higgsino and light 

Left and Right Sleptons 

LHLRS 

(Sec.4.2) 

u - Ml 
g 2 

mj = mj = (m~o + m~±)/2 

I'L IR A! Xi 

Light Higgsino and light 

Left Sleptons 

LHLS 

(Sec.5.2) 

u- Ml 

P 2 

mr = (m~o + m~±)/2,mr = 2 TeV 

‘.l v Xi Xi '' ir 

Light Higgsino and Heavy 

Sleptons 

LHHS 

(Sec. 6) 

2 

m~ = a+ 200 GeV 

l l,r 

xf-Mixed 

Light Mixed and light 

Left and Right Sleptons 

LMLRS 

(Sec.4.3) 

n = 1.05 M 2 

m~ = m~ = (m~o + m~±)/2 

Light Mixed and light 

Left Sleptons 

LMLS 

(Sec.5.3) 

n = 1.05M 2 

mr = (m ~0 + m~±)/2 1 mr = 2 TeV 

II Xi Xi Lr 


Table 1: Summaries of the models analyzed in this work. The parameter choice for each case is 
presented in the last column. For all the analyses we take M\ « M 2 , M\ ^ g to make the LSP 
predominantly a bino. Two representative values of tan /? = 6 and 30 are considered in this analysis. 
For the LHHS model, however, we consider only tan /3 = 30 case for reasons discussed in the text. 

DM relic density and (g — 2)^ in the higgsino models are also quite distinctive. All these points will 
be elaborated in the rest of this paper. 

We will further confront each APS thus obtained with other constraints like those from direct 
[37,38] and indirect [39,40] dark matter searches. We include current limits as well as compare 
our results in relation to future reaches of these experiments [39,41]. We also keep in mind the 
sizable uncertainties involved in these constraints (see Ref. [24] and the references therein). We 
emphasize that unlike the analysis of Ref. [24] , here we have a sizable amount of higgsino content 
within the LSP because g and M\ are not widely separated. This leads to a considerable increase in 
the spin-independent LSP-nucleon scattering cross-section. Additionally, a larger higgsino content 
within the LSP generically increases the spin-dependent LSP-nucleon cross-section. This, in turn, 
increases the gravitational capture cross-section of the LSPs within astrophysically dense regions like 
the core of the Sun. In addition, the LSPs with larger higgsino content may potentially contribute 
to neutrinos created within the Sun via LSP pair annihilation. Thus, the IceCube experiment puts 
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limits on the muon flux at the detector site which are not far away from the predictions of the 
models. 

We next consider an extended scenario in which only one strongly interacting sparticle is 
within the reach of the LHC-13/14 TeV experiments and assume this sparticle to be the gluino. 
The purpose is to study the feasibility of characterizing different higgsino models from their gluino 
decay signatures. We explore the gluino mass limits obtained at the LHC -8 TeV experiments via the 
n-leptons + m-jets (with or without b-tagging) + frp signals with different values of m and n [4-7]. 
By selecting a few benchmark points (BPs) with different characteristics we compute the revised 
gluino mass limits at the generator level by using PYTHIA (v6.428) [36]. This gives the sensitivity 
of various higgsino models to signals with different values of m and n and helps to anticipate the 
future search prospects. 

The plan of the paper is as follows. We note that because of the enhancement of higgsino 
components within X 3 and xf. the production cross-sections of xf-x® an d xf-Xz i n the xf- 
higgsino (xf-mixed) models are significantly (modestly) reduced from those in the xf-wino models. 
The characteristics of various cross-sections will be discussed in Sec. 2 with numerical examples. 
The relevant constraints for the higgsino and mixed models as well as the procedures for parameter 
space scanning and simulation of the trilepton signal will be discussed in Sec. 3. The parameter 
spaces allowed by the main constraints in the xf -higgsino and xf-mixed models characterized by 
different mass hierarchies among the sleptons and the elctroweakinos will be presented in Secs. 4-6. 
We will explore the prospects of direct and indirect detection of dark matter in Sec. 7 and Sec. 8 
respectively. The extended models with a light gluino and electroweak sparticles, introduced in the 
last paragraph, will be taken up in Sec. 9. Our main results and the conclusions are summarized 
in Sec. 10. 

2 Production of chargino-neutralino pairs in different models 

The size of the chargino - neutralino production cross-section with different neutralinos accompa¬ 
nying xf is significantly different in the wino, higgsino and mixed models. In the wino models the 
main signal, namely 3 l + fir comes from pp —>• Xixfi > since the production of the heavier neutralinos 
like X 3 or y® is highly suppressed. It may be recalled that xf and X 2 masses are degenerate in this 


case and are controlled by the soft breaking mass M 2 for the SU(2) gauginos while the masses of 
X 3 and Xa are governed by the higgsino mass parameter g, with fi » M 2 . This degeneracy holds 
for all m~± to a very good approximation. 


Masses 

PI 

P2 

P3 

and 

Model 

Model 

Model 

cross-sections 

Wino 

Mixed 

Higgsino 

Wino 

Mixed 

Higgsino 

Wino 

Mixed 

Higgsino 

m ~0 
x? 

150 

150 

150 

150 

150 

150 

150 

150 

150 

ra~± 

X? 

200 

200 

200 

300 

300 

300 

650 

650 

650 

m ~0 
x° 

201 

210 

219 

300 

304 

304 

650 

651 

652 

m ~0 

X3 

421 

269 

221 

604 

370 

312 

1256 

722 

657 

o-(PP-txfx°) 

0.7 

0.43 

0.194 

0.129 

0.083 

0.037 

0.00207 

0.00135 

6.9xl0 -4 

<T(PP -> xf X3) 

10" 3 

0.06 

0.209 

10~ 4 

0.011 

0.037 

- 

0.00022 

6.63xl0 -4 

&Total 

0.70 

0.49 

0.403 

0.129 

0.094 

0.074 

0.00207 

0.00157 

0.00135 


Table 2: Table showing the relevant masses and the cross-sections for three parameter points PI, 
P2 and P3. Here all the masses are in GeV and cross-sections are in pb. 

In contrast, the higgsino models are characterized by // << M 2 . As a result, X 2 an d X 3 are 
nearly mass-degenerate with xf where all masses are essentially determined by /r. The degeneracy 
is more exact as m~± increases. Here both X 1 X 2 an d X 1 X 3 production cross-sections are significant. 
In spite of this, the total chargino - neutralino production cross-section is smaller than that in a 
generic wino model with similar chargino and neutralino masses. In the mixed models (fj, ~ M 2 ), 
the cross-sections typically have intermediate values with respect to those of the wino and the 
higgsino models. The x^ and X 2 are nearly degenerate and the degree of degeneracy increases as 
m~± increases. There is always a much larger mass difference between X 2 an d X 3 compared to the 
higgsino model. It may be noted that these features are independent of Mi, the 1/(1) gaugino mass, 
as long as the LSP is bino-dominated, i.e, Mi « M 2 , //. The choice of tan /?, the ratio of the Higgs 
vacuum expectation values, only has a marginal impact. 

Production of x^ — X 2 /X 3 occurs through the process —>• xfx 2 3 ) where q and q' for 

the first two generations (i = 1,2) refer to up and down type of quarks respectively. When the 
first two generations of squarks are heavy, s-channel W boson exchange becomes the dominant 
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production process. We note that the — xf — X 2 3 coupling is contributed by terms involving 
products of wino components as well as terms having products of higgsino components of the 
relevant electroweakinos. The former terms typically dominate over the latter. Thus, as we move 
from the wino model to the higgsino model , the gradually diminishing wino contents of xf and X 23 
render smaller and smaller cross-sections for xf ~ X 2 /X 3 production. 

The above features are illustrated by Table 2. Here the LSP mass is fixed at 150 GeV. The 
parameters M 2 and /1 are differently chosen in the three models, in particular, /r = 2 M 2 (wino), 
H = M 2/2 ( higgsino ) and [i — M 2 (mixed). Throughout this paper we shall present the numerical 
results with these characteristic choices. It follows from Table 2 that for a wide range of the 
chargino mass the total chargino - neutralino production cross-section in the higgsino (mixed) 
model is 60-65% (70-75%) of that in the wino model. All next to leading order (NLO) cross-sections 
are calculated with PROSPINO 2.1 [42]. The reduction in the production cross-section is one of the 
reasons for relaxed mass limits in the higgsino and mixed models. However, the LHC mass limits 
are also sensitive to the choice of the slepton sector which will be addressed in Secs. 4, 5 and 6 . 

3 The Methodology 

In this section we summarize the constraints that we use to restrict the parameter spaces of several 
higgsino and mixed models. We also present brief sketches of the simulation using PYTHIA (v6.428) 
as well as the procedure for scanning the parameter space. 

3.1 The Constraints 

The three entries listed below are characterized by relatively small theoretical/experimental uncer¬ 
tainties and they constitute what we call the three major constraints. 

1. The LHC constraints from the chargino-neutralino searches in the trilepton channel used 
in this paper are from the ATLAS conference report [43] with C = 20 fb -1 data, which is 
the source of the published paper [16]. Similarly, for the constraints in the slepton sector 
we use the conference report [44], It may be noted that there is no major difference in 
the exclusion contours among the published versions and their predecessors. The ATLAS 
collaboration quotes the upper limits on the number of events in any new physics model 
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at 95% CL subject to different sets of selection criteria (see the next subsection). We have 
simulated the trilepton signal using the same sets of kinematical selections at the generator 
level using PYTHIA (v6.428) for the higgsino and mixed models analysed in this paper. The 
model independent limits then enable us to sketch the exclusion contours in the m~± — m ~o 
plane in different models for representative choices of the other parameters as detailed in 
the following sections. As already noted in the introduction the slepton constraints hardly 
depend on the composition of the charginos and neutralinos heavier than the sleptons. We 
have, therefore, directly used the exclusion contours from light L and LR slepton searches in 
the corresponding wino models obtained in [24] . 

2. The anomalous magnetic moment of the muon (a^ = \{g — 2)^) [45] is an important probe for 
new physics beyond the standard model (SM). There is a significantly large deviation (more 
than 3a) of the SM prediction [46,47] from the experimental data [35]. Contributions to 
can be categorized into three parts : a part coming from pure quantum electrodynamics, 
electroweak contributions and a hadronic part. SUSY contributions to a^, namely a® USY , 
scale with tan /3. It can also be large when chargino, sneutrino, neutralino and srnuons are 
light [48]. Thus, it is possible to constrain the SUSY parameter space effectively with given 
upper and lower limits of A= a® xp — a® M . The deviation of the experimental data from 
the SM calculation amounts to [47] 

Aa M = a“ p - = (29.3 ± 9.0) x 10~ 10 . (1) 

A partial list of analyses regarding a® USY in SUSY models is provided by Refs. [48-51]. With 
the Higgs mass at 125 GeV and stringent lower bounds on squark-gluino masses coming from 
the LHC, simplified models like rnSUGRA have become rather inefficient to accommodate 
the (g — 2)^ anomaly [52]. However, non-universal models can still successfully explain the 
above range of A [53]. It should be noted that the (g — 2)^ constraint is able to impose 
definite upper and lower bounds on the sparticle masses [54] . 

3. Following Ref. [24], relic density limits from the WMAP/PLANCK [29,30] is taken as, 

0.092 < Q^h 2 < 0.138. (2) 
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Apart from the above direct constraints we will further analyze our results in relation to the 
following dark matter detection limits. 

• The direct detection bound on spin-independent (SI) LSP-proton scattering cross-section a^ p 
is imposed using the LUX [38] data. 

• DM indirect detection constraints like the bounds on spin-dependent (SD) LSP-proton scat¬ 
tering cross-section cr-^ and limits on muon flux given by the IceCube [39] are also important 
in our case. 

We use SuSpect version 2.41 [55] for spectra generation and a® USY calculation. SUSYHIT [56] is 
used for obtaining the decay BRs of the sparticles. DM relic density and observables related to its 
direct and indirect detection are computed using micrOMEGAs version 3.2 [57]. 

3.2 The Simulation 

In Ref. [43] the ATLAS collaboration defined six signal regions (SRs) : SRnoZa, SRnoZb, SRnoZc, 
SRZa, SRZb and SRZc. Table 1 of Ref. [43] includes the details of the cuts corresponding to each 
signal region. The corresponding upper limits on the number of new physics events are listed in 
Table 4 of Ref. [43]. In all models analysed in this paper we simulate the trilepton signal for the 
above SRs with a given m~± by increasing the LSP mass in small steps. Below a certain LSP mass 
the above ATLAS upperbounds on the number of new physics events are violated for at least one 
SR. A point on the exclusion contour is determined in this way. On the other hand when m~± 
is varied, all the LSP masses are allowed above a certain m~± which is the lower limit on this 

Xi 

parameter. The validation of our simulation and other details have been described in Ref. [24] . We 
use the same setup for the present work. 

3.3 Scanning the parameter space 

We take the strong sector parameters to be heavy by choosing M 3 = 2 TeV and fixing the masses of 
the first two generations of squarks and the mass (Ma) of the pseudoscalar Higgs boson at 3 TeV. 
The mass parameters for the third generation of squarks are fixed at 1.2 TeV. The top trilinear 
parameter At is varied in the range -5 TeV < A t < 5 TeV in order to obtain the lighter Higgs scalar 
mass to be in the interval 122 < < 128 GeV. All the other trilinear parameters except At are 
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taken to be zero. A theoretical uncertainty of 3 GeV in the computation of Higgs mass is considered 
here. This spread in arises from the uncertainties in the higher order loop corrections up to 
three loops, that due to the top-quark mass, renormalisation scheme, scale dependence etc [58]. The 
Higgs bosons other than the lightest one are assumed to be decoupled. We perform scan over the 
parameters Mi, and M 2 in the range 100 GeV-1 TeV. The slepton masses are correlated to m ~0 and 
m~±/m ~0 in this work and the nature of correlation are described in the appropriate subsections. 
Values of the relevant SM parameters are taken as m^ ole = 173.2 GeV, m^ IS = 4.19 GeV and 
m T = 1.77 GeV. Finally, we consider only the positive sign of /a in this analysis. 

4 The xf _ Higgsino and xf~Mixed Models with Light Left and 
Right Sleptons 

In this sequel analysis that extends our previous work [24] on xf~wino models towards x^-higgsino 
ones we will review a few salient points as well as refer to some figures of the earlier analysis for the 
sake of clarity and easier understanding of the present work. The analysis of various x^-higgsino 
and x^-mixed models using the methodology sketched in the last section, will be presented in this 
and the next two sections. 

4.1 A brief review of the Light Wino and light Left and Right Sleptons (LWLRS) 
model 

In the Light Wino and light Left and Right Sleptons (LWLRS) model [24] , the L and the R types 
of sleptons were assumed to be mass degenerate (modulo the D-term contributions) with a common 
mass : (aq M\ + aq M 2 ). Three choices were considered: i) aq = aq = 0.5, ii) aq = 0.25,aq = 0.75 
and iii) aq = 0.75, aq = 0.25, with the slepton mass lying between and M 2 . The models with 
aq 7 ^ aq are referred to as tilted models denoted by the LWLRS~±(LWLRS~o) with the slepton 
mass closer to the xt (xi) mass. The ATLAS group did not interpret their data in any form 
of the LWLRS model. However, as we shall show in the next subsection, some versions of the 
corresponding higgsino model are indeed rather intriguing since they are practically unconstrained 
by the LHC data. 

The impact of the three major constraints (see Secs.l and 3) on the parameter space of the 
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LWLRS model for low tan/? (= 6) may be seen in Figure 4a of Ref. [24]. For a negligible LSP 
mass, the acceptable value of m~± is above 610 GeV. On the other hand, for higher values of the 
LSP masses the limits that had been obtained primarily from the trilepton searches earlier, became 
weaker. There were two branches in the parameter space consistent with the WMAP/PLANCK 
data [29,30]. This feature is common to most of the wino models with low tan/3. The fact that 
relic density limit is satisfied in the upper branch is related to LSP-sneutrino coannihilation. The 
lower branch is consistent with the relic density limits via LSP pair-annihilation into the Higgs 
resonance, corresponding to m~ o ~ m^/2 9 . In spite of satisfying the DM relic density limits, a 
large part of the above parameter space corresponding to low m~± is disfavoured by the trilepton 
and slepton searches. This indeed is the case in all the wino models with x\ = X 2 = 0.5. However, 
the constraints from slepton searches are significantly relaxed in the LWLRS~± model and one 
obtains m~± > 450 GeV for negligible LSP masses(see Fig. 6a of Ref. [24]). The results of Ref. [24] 
derived from slepton searches are being readily adapted in the corresponding higgsino or mixed 
model of the present analysis for reasons already discussed in Sec. 1. However, the entire parameter 
space allowed even by the relaxed constraints from slepton searches is consistent with the (g — 2)^ 
data at best at the 3er level. This tension at low tan/3 from the (g — 2)^ constraint exists in all the 
wino models as already been noted [24], This feature is also shared by the higgsino models as we 
shall see in the following sections. 

A choice of a large tan /3 (= 30) in the LWLRS model hardly changes the LHC constraints, as 
may be seen in Fig. 4b of Ref. [24]. However, a significant amount of parameter space is discarded 
because f\ becomes the LSP as a result of enhanced degree of Left-Right mixing considering a large 
value of tan /3. On the other hand, since f\ is significantly lighter than the other sleptons, r-rich 
final states deplete a part of the trilepton signal 10 . This depletion, however, is more effective in 
the parameter space beyond the reach of the LHC Run-I experiments. We should here point out 
that the small mass difference between f\ and Xi which results into satisfying DM relic density via 
LSP-stau coannihilation (see the upper branch of the region allowed by the WMAP/PLANCK data 
in Fig. 4b of Ref. [24]) makes the tau-tagged signatures difficult to observe. The Higgs resonance 

9 In some of the cases the Higgs resonance region is accompanied by a small Z resonance annihilation branch. 

However, we will mostly focus on the former in our present analysis. 

10 The potential of final states with multiple rs from x/ an d X 2 decays as new signals at the LHC was emphasized 

in Ref. [59] . 
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region on the other hand disappears due to a vanishingly small higgsino component of Xi which 
severely suppresses the hxiXi coupling at high tan/3, a generic feature shared by all the wino 
models with a bino-like LSP. As expected the tension with the (g — 2 )^ constraint is relaxed in this 
high tan /3 scenario. As a result there is a narrow APS consistent with all the major constraints. 

4.2 The Light Higgsino and light Left and Right Sleptons (LHLRS) model 

In this subsection we focus on the models characterized by Light Higgsino and light Left and 
Right Sleptons (LHLRS), where the sleptons are assumed to have masses nearly halfway between 
the masses of the lightest neutralino and the lighter chargino 11 while xf and X 2 are higgsino- 
dominated in nature. In this scenario with Mi < u < M 2 , m ~±, m~ 0 and m~o are close to u. In 

Xl X2 X3 

this paper we restrict ourselves to the regions of parameter space where the LSP is either bino- 
dominated (Mi << g) ora strong admixture of bino and higgsinos with the former field dominating 
over the others. In the latter case, the LSP and xf are almost mass degenerate. 

We do not however include the scenario where the LSP is purely a higgsino or higgsino- 
dominated (/r ^ Mi). In this parameter space all the lighter electroweakinos including the LSP 
are higgsino-like and nearly mass degenerate. As a result, no interesting LHC signature, except 
for the well-known monojet + pr signal, is expected. DM searches by the LHC collaborations in 
this channel yielded only null results [60]. However, the results have not been interpreted in terms 
of the pMSSM scenarios. According to the analysis of Ref. [61] using the LHC Run-I data the 
bound on the LSP mass in this case is rather weak. Moreover, for the electroweakino mass ranges 
of current phenomenological importance, the DM is found to be underabundant when g, ^ Mi [62] . 

In Fig. 1(a), we present the results of the LHLRS analysis for tan/3 = 6. The brown and 
the green regions signify the parameter space consistent with the (g — 2) M constraint at 3a and 
2 <t levels respectively. The magenta line represents the reference contour which is the strongest 
m ~0 — m~± mass limit as obtained in the Light Wino and light Left Slepton (LWLS) model [see 
Fig. la of Ref. [24]]. The black line is the contour obtained by our simulation for the LHLRS model 

11 Placement of slepton masses follow similar relationship with the masses of the LSP and the lighter chargino as 
mentioned in the previous subsection. This is also accompanied by similar tilted scenarios like the LHLRS~-t and 

Xl 

LHLRS—q models where sleptons are closer in masses to that of xf and Xi respectively by specific amounts. 

Xi 
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(a) (b) 

Figure 1: The brown, green and yellow regions show the parameter space consistent with the 

(g — 2)^ constraint at 3a, 2a and lcr level respectively. The magenta line is the reference contour 
which represents the strongest m ~o — m~± mass limit as obtained in the Light Wino and light Left 
Slepton (LWLS) model [see Fig. la of Ref. [24]]. The black line is the contour obtained by our 
simulation using ATLAS data. The black dashed line represents the exclusion contour from ATLAS 
slepton searches. The grey region to the left and above the coloured portion of the parameter space 
is either theoretically discarded or disallowed by our requirement of having a bino-dominated LSP 
or disallowed for the requirement of Xi to be the LSP. The red points satisfy the WMAP/PLANCK 
data on DM relic density. The m~± = m ~o line is also shown for the plot in the left panel. For the 
plot in the right panel, this line does not exist since f\ becomes the LSP for m~± ~ m ~o because 

Xi X\ 

of large mixing. 

using ATLAS data. The black dashed line represents the exclusion contour from ATLAS slepton 
searches. The grey region to the left of the APS is theoretically disallowed. The small change in 
this region compared to that of Fig. 4a of Ref. [24] is due to the choice of a smaller stop mass 
(here mq: — 1 TeV, the same in the previous analysis was 2 TeV) and its effect on the electroweak 
symmetry breaking (EWSB) scale, Q = m~ , as used in the code SuSpect [55]. The grey 
region just below the m~± = m~ o line corresponds to /j, ^ M\ which leads to a higgsino-dominated 
LSP, a scenario kept out of the domain of our analysis for reasons already discussed. The grey 
region above the aforesaid line is discarded because Xi is chosen to be the LSP, a candidate for 
DM. We should mention that the theoretically disallowed regions which we will often quote in this 
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work arise from the facts mentioned above. We also take into account the LEP limits [63] for 
sparticle masses, if necessary. 

It is apparent that apart from a small region of the parameter space corresponding to low m ~±, 
the collider limits are relaxed significantly with respect to the wino dominated scenario (Fig. 4a 
of Ref. [24]). For negligible m~ o, a substantially weaker bound (m~± 380 GeV) compared to 

the wino model is obtained via the trilepton signal. However, the slepton mass limits discard m~± 
up to 600 GeV. On the other hand, with a modestly increased value of m~ o such as 150 GeV, 
the lower bound on m~± is about 330 GeV in order to satisfy both the limits from trilepton and 
slepton searches. For larger LSP masses, the chargino mass bounds are even weaker, as expected, 
and eventually disappears for m~ o ^ 200 GeV. In addition to the suppression of the total chargino- 
neutralino production cross-section as discussed in the Sec. 2, the decrease in the leptonic BRs of 
the electroweakinos reduces the trilepton signal. For small m~±, the decays like xf W±Xi and 
X 2 /X 3 —> h/Z Xi are kinematically forbidden. Hence, the number of two-body leptonic decays of 
xf and X 2 remain practically unaltered compared to the wino model. Moreover, the lepton pairs 
from the decays of xfl compensate the overall decrease in the cross-section and this leads to almost 
unchanged LHC limits. The situation changes significantly when the decays xf W 1 * 1 x? an d 
X 2 / X?i —>• h/Zxi are kinematically allowed. The BR(x^ —>• W± Xi) being dominant (~ 75%), 
the overall lepton fraction drops down as W-boson decays with the usual small leptonic BRs 12 . 
Moreover, x .2 an d X 3 decay dominantly into the gauge/Higgs boson channels, thus depleting the 
trilepton signal further 13 . 

Regarding the constraint from (g — 2)^, the major contribution comes from the neutralino- 
smuon loop processes. However, the contribution from the chargino-sneutrino loop diagrams are 
hardly ignorable. In this low tan/3 scenario of Fig. 1(a), the (g — 2)^ constraint which typically 
is satisfied only at the level of 3<r. As expected, the SUSY contribution to (g — 2) M falls on the 
lower side of the 3 a limit rather than on the higher side. In other words, keeping in mind the 
level of uncertainties of both the SM prediction and the experimental data [45-47] , the above 3 a 
zone is indeed close to the zone where agreement with the experimental data is at higher level of 

12 This happens in spite of the presence of light sleptons which do not couple favourably to the higgsino like xf 
and x° 2 - 

13 Numerical values of the relevant BRs for selected BPs will be provided in Sec. 9. 
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confidence. The situation would change with more accurate SM predictions and more precision in 
the experimental front. We simply like to comment here that accepting such a relaxed level of the 
(g — 2)^ constraint may be worthwhile if there is a scope to open up a region of parameter space 
having important LHC signature(s). As discussed above, in this 3 a allowed parameter space, the 
decays xt Wx? and X 2 ^Xi would potentially lead to novel collider signatures like Whxi 
during the LHC Run-II. 

There are two separate branches consistent with the WMAP/PLANCK data in Fig. 1(a). In 

the parameter space with small slepton masses, bulk annihilation 11 of x^ may partly serve as the 

cause to satisfy the relic density limits via enhanced LSP pair-annihilation. But, this region of 

the parameter space is disfavoured by the LHC constraints anyway. In the upper branch, g is 

close to Mi so that the bino-like LSP has a sizable higgsino component. Thus, here one finds x]* 1 

mediated LSP pair annihilations to W + W~ to play a significant role. Annihilation into ZZ, Zh and 

tt through virtual Z exchange open up for higher LSP masses and also contribute significantly. In 

contrast, the mechanism for satisfying relic density limits in the similar branch of the wino model 

is mainly due to coannihilations of the LSP with stau/sneutrino (see Fig. 4a of Ref. [24]). It is 

interesting to note that the DM constraint also provides with bounds on the sparticle masses from 

above since the upper branch ends abruptly above a certain value of m~ 0 . This is due to the fact 

that for large m~±, the chargino mediated LSP pair annihilation cross-section is suppressed. In 
Xi 

addition, the annihilation cross-section to tt pairs goes as (^-) 2 [26,64], which becomes inefficient 

x ? 

due to increasing value of m~ 0 , resulting into over-abundant DM. 

The DM relic density satisfied lower branch forming almost a line parallel to the m~± axis 
corresponds to LSP pair-annihilation via s-channel Z/h resonance and it mostly falls on the 3u 
zone of the (g — 2) M constraint as discussed above. Moreover, a large portion of this line is excluded 
by slepton searches for relatively low m ~±. The constraints from slepton searches are considerably 
relaxed in the LHLRS-xf model as discussed in Sec. 4.2. 

Thus, novel LHC signatures via the decay x.2 hXi is a potential discovery channel at the 
LHC (see Sec. 9) provided we accept the 3a level of agreement with the (g — 2) M data. 

In Fig. 1(b) we present the results for tan /3 = 30. A larger value of tan /3 causes the lighter stau 
(fi) to become the next to the lightest supersymmetric particle (NLSP) due to larger L-R mixing. 

14 LSP-pair annihilation that occurs via t-channel slepton exchange. 
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We note that in comparison to a similar wino-dominated xf scenario like what appears in Fig. 4b 
of Ref. [24], the theoretically discarded region corresponding to low m ~o is smaller in Fig. 1(b). It 
is partly due to the reasons already discussed for Fig. 1(a). Moreover, /r being relatively low, as 
demanded by a higgsino dominated xf, the mixing in the stau sector is suppressed. As a result, 
the portion of the parameter space with stau as the LSP is smaller. No m~± = m ~o line exists in 
this case as f\ becomes the LSP for m~± ~ m~o because of large mixing. Fig. 1(b) shows further 
that the LHC exclusion contour from the trilepton data is too small to be of any importance in 
determining the APS. In fact, the LHC limits are fully encompassed by the theoretically excluded 
region. In addition to the reasons already discussed, the constraints are also weakened because of a 
significantly reduced mass of f\ in comparison to the same of sleptons of the first two generations. 
Consequently, there is a significant suppression in the BRs of the decays of xf, X 2 an d X 3 into 
states involving e and/or fj.. This explains the shrinkage of the LHC forbidden region compared to 
Fig. 1(a). 

The (g — 2)^ constraint that typically scales with tan /3 is satisfied much easily in this large 
tan f3 scenario even at 1 a level. The largest contribution to (g — 2)^ is from the neutralino-smuon 
loop. But the chargino-sneutrino loop also has a significant contribution. 

The features of the upper branch allowed by the WMAP/PLANCK data are similar to those in 
Fig. 1(a). However, since f\ is the NLSP in this case stau coannihilation also becomes important for 
higher LSP masses. Nevertheless, the upper branch is also truncated in this case, thereby imposing 
an upper limit on the sparticle spectra. This bounded region is consistent with both the (g — 2)^ 
and the LHC constraints. We recall that in the wino dominated scenario (see Fig. 4b of Ref. [24]), 
the Higgs resonance region disappears at large tan /3 due to vanishingly small higgsino component 
of Xi which severely suppresses the hxfxi coupling. In contrast, /r being smaller in the present 
scenario, xf though dominantly a bino, has sufficient higgsino components. Thus, the h-resonance 
region extends to large values of m~±. The bulk of this region, though allowed by the (g — 2)^ 
constraint, is disfavoured by the direct slepton search data. However, as already noted, this region 
opens up in the LHLRS~± model. 

Xi 
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(a) (b) 

Figure 2: Plots in the m~± — m ~o plane for the LMLRS scenario with tan B=6 and 30. Colours 

Xi Xi 

and conventions are the same as in Fig. 1. 


4.3 The Light Mixed xt with light Left and Right Sleptons (LMLRS) model 

The analysis of this subsection is based on the Light Mixed xf with light Left and Right Sleptons 
(LMLRS) scenario where the states xf, X 2 an d X 3 are admixtures of large wino and higgsino 
components while the slepton masses are chosen as in the previous subsection. We consider p = 
1.05 M 2 , a choice that corresponds to a large higgsino-wino mixing while keeping xf, X 2 an d X 3 
to be still dominated by the higgsinos. On the other hand, x.i ma y have a significant amount of 
higgsino component in some region of parameter space, although it is by and large a bino-dominated 
state (/r > Mi). We analyze the production processes : pp —>• xf X^/xls- In Fig. 2(a), we show our 
results for tan (3 = 6. 

A gradual weakening of the LHC constraint from the trilepton search is evident from the 
LWLRS model (Fig. 4a of Ref. [24]), to the LMLRS model (Fig. 2(a)) and then to the LHLRS 
model (Fig. 1(a)) while the higgsino fraction within xf increases steadily. This is as expected from 
the discussions on the cross-sections in Sec. 2. 

Thus, for a negligible LSP mass the lighter chargino mass limit that is allowed via the trilepton 
search is about 550 GeV in Fig. 2(a) compared to about 390 GeV in Fig. 1(a). However, the model 
independent results of slepton searches push the above limits in the figures to almost an identical 
value (~600 GeV). 
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34 
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Table 3: The sparticle spectra corresponding to different points chosen from Fig. 2 and Fig. 4. All 
the masses are in GeV. 

We observe that in the small m~± region of Fig. 2(a), the limits from the LHC trilepton searches 
are practically unaltered with respect to the corresponding wino model (Fig. 4a of Ref. [24]), in 
spite of differing wino and higgsino components within xt between the two scenarios. Here BR 
of X 2 decaying into invisible vv final state reduces significantly due to an increase in its higgsino 
fraction. This on the other hand is compensated the reduction in the cross-section via an increase 
in the leptonic BR leading to an almost unchanged LHC limit. We further note that X 3 does 
not contribute to the trilepton signal appreciably since it decays principally via Zx\ mode with 
BR around 80%. Moreover, the production cross-section for X3 xt is rather small compared to 
X2 xt since m~ 0 is somewhat larger than m~o. Some of the above features are illustrated by two 
representative points taken from Fig. 2(a)) in Tables 3 and 4. In the region of high m~±, the decay 
modes x .2 ~^Xi and xf ► W^xt open up (the BRs lie around 10%). X 3 decays dominantly into 
Zx 1 as before but a sizable fraction goes to hxt Thus, the LHC exclusion contour shrinks when 
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Decay Modes 
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Table 4: The decay modes and BRs of different points taken from Fig. 2 and Fig. 4. 

compared with Fig. 4a of Ref. [24]. These features are illustrated in Table 4. 

Regarding (g — 2) M , in all the cases the LHC allowed parameter space is consistent at 3cr level. 
Both the loops involving neutralinos as well as charginos contribute with comparable magnitudes. 

Along the upper red dotted line, the main contribution to the observed relic density are chargino 
mediated LSP annihilation into the final state W + W~ and LSP pair annihilation into tt . Some 
amount of bulk annihilation is also present for low m~± although this region of parameter space is 

Al 

discarded by the LHC constraints. Unlike the higgsino model of Fig. 1(a) the upper branch does 
not end abruptly (at m~± & 400 GeV). This is due to the fact that Xi ~ xf ~ coupling which 
is behind the xf mediated annihilation to W + W~ is stronger because of larger wino content of the 
LSP in the mixed model. In the corresponding wino model (Fig. 4a of Ref. [24]), however, the relic 
density falls within the observed range mainly due to stau/sneutrino coannihilations. The lower 
branch of red points consistent with the observed DM relic density arises through Z/h resonance. 
Fig. 2(b) illustrates the results for tan/3 = 30. The small change in the theoretically discarded 
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region is due to reasons explained earlier (see Sec. 4.1). Here the parameter space excluded by the 
LHC constraint shrinks with respect to that in the wino model (Fig. 4b of Ref. [24]). The reasons 
are the same as the ones discussed in the context of Fig. 2(a) and 1(b). However, X 3 decay has no 
bearing on the LHC limits in this case. Significantly larger parameter spaces consistent with the 
(g — 2)^ constraint at the level of la and 2 a are available in this high tan /3 scenario. On the other 
hand, for the DM constraint as satisfied in the the upper red dotted branch, the main DM producing 
mechanisms are stau-LSP coannihilation and stau-stau annihilation as in the wino model (Fig. 4b 
of Ref. [24] ). Chargino-mediated annihilation and annihilation into tt pairs are also present albeit 
to a lesser extent. The features of one representative point in this parameter space is illustrated in 
Tables 3 and 4. Here the Higgs resonance strip extends to higher m~± compared to the wino model 
(see Fig. 4b of Ref. [24]) where it was practically absent. This is due to a modest increase in the 
higgsino component of the LSP. However, this parameter space is strongly disfavoured by the LHC 
limits even in the LMLRSx ^ 1 model. 

We end this section by noting that in both the LHLRS and LMLRS models with high tan /3 
we obtain several APSs consistent with the main constraints. 

5 The Higgsino and Mixed Models with Light Left Sleptons 

5.1 A brief review of the Light Wino and light Left Sleptons (LWLS) model 

We recall that in Fig. 1 of our previous analysis (Ref. [24]), the model characterized by Light 
Wino and light Left Sleptons (LWLS) yields the strongest mass limit on xt- For negligible m~ 0 , 

A-l 

m~± ranging up to 610 GeV is excluded. In all cases the low tan (5 scenarios are consistent with 
the (g — 2) m constraint at best at the 3er level. On the other hand, for large tan (5, the SUSY 
prediction has a better agreement with constraint which is satisfied at la or 2 a levels. In all 
cases the parameter space consistent with the WMAP/PLANCK data has two limbs. In the upper 
limb various coannihilations lead to the desired relic density. The lower limb which represents the 
LSP pair annihilation through the h-resonance is either absent or strongly disfavoured by the LHC 
constraints in the high tan (3 scenarios. 
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5.2 The Light Higgsino and light Left Sleptons (LHLS) model 


In this section we replace wino of the previous subsection by a higgsino and focus on models with 
Light Higgsino and light Left Sleptons (LHLS) where xt an d X 2 are higgsino-dominated and left 
slepton masses are kept midway between m~± and m~ 0 . Right sleptons are taken to be heavy. 

Fig. 3(a) shows our results for the model with tan [3= 6. If we compare this with Fig. la of 
Ref. [24] we see that the theoretically discarded regions are almost the same in both cases. The 
small change is due to a different choice of EWSB scale as noted earlier. The LHC limits from the 
trilepton searches are significantly degraded in the higgsino case. The reasons are the same as the 
ones in Sec. 4.2. For low LSP masses, however, one finds a region with m~± 600 GeV to be 

allowed from slepton search. This limit will be relaxed to 410 GeV in the tilted LHLS~± models 
(See [24] Figs. 3a and 3b). 

The major contribution to (g — 2)^ comes from the chargino-sneutrino loop as expected in 
models with light L-type sleptons (see Ref. [24] for details). The (g — 2)^ allowed regions at the 
level of lcr, 2 a and 3<r are similar to those of Fig. la of Ref. [24]. The observations regarding the 
points satisfying the WMAP/PLANCK limits are similar to those of Fig. 3(a) i.e. the lower red 
dotted branch comes from Z/h resonance annihilation and the upper branch from annihilation to 
W + W ~, tt, ZZ pairs. There is a small amount of sneutrino coannihilation in the upper branch 
along with a little bulk annihilation for very low m~± disfavoured by the LHC bounds. The reasons 
for the abrupt end of the upper branch indicating upper bounds on m~± and m~ 0 are already 

Xi Xi 

discussed in Sec. 4.2. The (g — 2)^ constraint is satisfied only at the level of 3<r. 

Our results for tan f3 = 30 is shown in Fig. 3(b). The collider exclusion limit is the same as that 
of Fig. 3(a). Since the right sleptons are heavy, there is no large mixing in the stau sector. As a 
result, the LHC exclusion contour remains unaltered in spite of the change in tan j3. The dominant 
contribution to (g — 2)^ comes from the chargino-sneutrino loop diagram and it is characteristically 
similar to Fig. 3(a) except that there is the usual tan/3 enhancement leading to valid lcr and 
2 <t regions. Annihilation and coannihilation mechanisms of the LSP in the regions satisfying the 
WMAP/PLANCK constraints in this scenario are similar to those of Fig. 3(a). The LSP pair 
annihilation into the Higgs resonance is allowed in the LHLSx]* 1 model. 
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(a) (b) 

Figure 3: Plots for the LHLS model with tan j5= 6 (a) and 30 (b). Colours and conventions are 

the same as in Fig. 1(a). 


5.3 The Light Mixed xf and light Left Sleptons (LMLS) model 

In Fig. 4(a) we show our results for the Light Mixed xf and light Left Sleptons model (LMLS) with 
tan /3= 6. The collider exclusion contour is weakened compared to the reference contour of Fig. la 
of Ref. [24], The reasons are almost same as those mentioned in the discussions for Fig. 2(a). Major 
contribution to (g — 2)^ is provided by the chargino-sneutrino loop diagram. The points satisfying 
the WMAP/PLANCK limits in the lower branch arise as a result of Z/h resonance annihilation 
processes. For the upper branch, main mechanisms are annihilations into W + W ~, ZZ , tt pairs. 

Fig. 4(b) shows the results for large tan /3 in the LMLS scenario. The observations regarding 
the upper limb consistent with the WMAP/PLANCK constraints are similar to those of Fig. 4(a). 
As in all high tan f3 scenarios there are portions of the parameter space consistent with the (g — 2)^ 
constraint at la and 2a levels leading to an APS consistent with all the major constraints. We 
show the features of one representative point for this scenario in Tables 3 and 4. 
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(a) (b) 

Figure 4: Plots for the LMLS model with tan/3=6 (a) and 30 (b). Colours and conventions are the 

same as in Fig. 1(a). 



Figure 5: Plot for the LHHS scenario with tan/3 =30. The magenta line represents the exclusion 
contour at 8 TeV by the ATLAS collaboration [43] for the wino model. The black line stands for the 
much weaker exclusion contour in the corresponding higgsino model obtained by our simulations. 
Other colours and conventions are the same as in Fig. 1(a). 
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6 The Light Higgsino and Heavy Sleptons (LHHS) model 


In the Light Higgsino and Heavy Sleptons (LHHS) scenario the chargino decays to Wxi with 
100% BR and the decay X 2 ^Xi has 100% BR for ( m~o — m~ 0 ) < rrih■ For larger m~o, the mode 
X 2 hx°i opens up and its BR can be « 40% near threshold. With increasing value of m~o — m~ 0 , 
the BR (%2 —>• hxi) may be as large as 60%. The presence of this mode is the main difference with 
the LWHS scenario and leads to the (W + h + ffip) signal. The BR of the mode x .3 ^Xi is around 
70%. However, the LHC limits in this case are degraded due to a reduction in cross-section as 
already discussed. For m~± > 200 GeV, all the LSP masses are allowed. On the other hand, for 
negligible LSP masses, a weaker limit (m~± ^ 175 GeV) is obtained. 

The above features, and, consequently the LHC exclusion contours, are fairly independent of 
tan /?. However, as already shown in the last two sections, the constraints from (g — 2)^ is effective 
for high tan /3. This is particularly true in a situation with heavy slepton masses that potentially 
reduces a® us ^ . Thus, only the case of tan (3 = 30 is shown in Fig. 5. The SUSY contribution to 
(g — 2)^ is dominated by the chargino-sneutrino loop diagram. The upper branch of the region 
consistent with the WMAP/PLANCK data arises due to chargino-mediated LSP pair annihilation 
to W pairs. There is also annihilation to fermion anti-fermion pairs through virtual Z exchange 
and some amount of LSP-chargino coannihilation. The Z/h-resonance region is also allowed in this 
scenario by the major constraints and the (g — 2)^ constraint is satisfied even at the level of la. 

An analogous case of a Light Mixed and Heavy Sleptons (LMHS) scenario would have the 
corresponding LHC exclusion contour lying in between the magenta and the black lines of Fig. 5. 
This is expected from the chargino-neutralino production cross-sections as discussed in Sec. 2. We 
do not present the details here since no qualitatively new feature of the signal emerges from this 
analysis. 

7 Direct Detection via Spin-independent scattering 

Spin-independent (SI) interaction of the lightest neutralino with quarks inside the detector nucleus 
occurs via s-channel squark exchange and t-channel Higgs exchange processes. With the present 
LHC bounds squarks are considerably heavy. Hence, the Higgs exchange diagrams would domi¬ 
nantly contribute towards 0 ^, the spin-independent x~ V scattering cross-section [65]. In this 
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context we note that the h(H)x 1 X 1 coupling involves products of the gaugino and the higgsino 
components of the neutralino diagonalizing matrix [3]. Unless \M\ — |/i|| & Mz, in the Higgs de¬ 
coupling zone [3] satisfying run — tua >> Mz, the Higgs couplings to a bino-dominated LSP is 
approximately given as below [66]. 


m z swtw r „, . . 

C hxx - M 2_ m 2 [Mi + nsm2p] : 

n mzswtw r, a 

°HXX ~ iT72- —V COS 2(3. 


( 3 ) 


M'f - fj? 

Here sw = sin Ow etc. with 9w as the Weinberg angle. Similar results for a wino or a higgsino- 
dominated LSP may be seen in Ref. [66]. Clearly, the above shows that the SI cross-section <T^ 
would be large when there is a significant amount of bino-higgsino mixing i.e. M\ ~ /x. This is unlike 
a pure gaugino or a pure higgsino DM when the associated SI cross-section would become quite 
small. We note that in contrast to our previous analysis [24] , where we considered a bino dominated 
LSP with n » Mi, M 2 , the present work has a significant amount of higgsino component in the 
LSP. As a result, is typically larger than what was seen in Ref. [24] and in most of the cases 
its values lie above the LUX [38] limit. However, we must note that there still exists a significant 
amount of uncertainty in the theoretical estimate of the SI cross-section (for a brief discussion see 
Ref. [24] and references quoted therein). This at least relates to issues like uncertainties in the 
determination of the strangeness content of nucleon, local DM density and velocity distribution 
profiles. All these uncertainties may accommodate lowering of the cross-section by an order of 
magnitude. We will present our results in the following subsections for the higgsino and the mixed 
models. We like to point out that no tilted scenarios have been included in our analyses on direct 
and indirect detection of dark matter. 


7.1 LHLRS and LMLRS 

The results for the SI direct detection are shown in Fig. 6(a) and Fig. 6(b) corresponding to i) 
higgsino-dominated xf (LHLRS) and ii) wino-higgsino mixed Xa (LMLRS) analyses of Fig. 1 and 
Fig. 2 respectively. Fig. 6(a) combines the results in the m ~0 — a ^ plane corresponding to two 
values of tan {3 (blue and cyan for tan/3 = 6 and 30 respectively) as used in Fig. 1(a) and Fig. 1(b). 
Similarly, Fig. 6(b) details the results corresponding to Fig. 2(a) and Fig. 2(b) for the same values 
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of tan /3 as mentioned above. Figs. 6(a) shows only the allowed points that satisfy the relic density 
limits, the (g — 2)^ constraint (upto 2a for large tan /3, and 3 a for small tan f3 cases) and the collider 
limits for m~o between 200 and 350 GeV. a^ for the points exceed the LUX limit, while staying 
within an order of magnitude of the same limit 15 . In regard to the mixed model (Fig. 6(b)), we 
similarly find that parameter zones satisfying m~o > 200 GeV are consistent with the major con- 
straints. Here also the points exceed the LUX limit by a similar amount. This is unlike the analysis 
of the LWLRS scenario of Ref. [24] where it was not difficult to satisfy the LUX limit. Since the 
parameter points correspond to a deviation below an order of magnitude from the LUX limit, we 
consider them to be presently acceptable in view of the uncertainties discussed before. There are 
a few points in the Higgs resonance region for the low tan /3 case. Since Higgs-pole annihilation 
occurs for m~ o ~ and the allowed points from this region correspond to very high m~± (see 
Fig. 2(a)), the LSP in this region is highly bino-dominated, thus having very small values of a^. 
The XENON1T experiment will conclusively probe these models. 




(a) (b) 

Figure 6: Direct detection results for the LHLRS (a) and LMLRS (b) scenarios with tan/3 = 6 

(blue) and 30 (cyan). The LUX and XENON1T limits are shown as black and red lines respectively. 
The points satisfying the WMAP/PLANCK and LHC limits are shown. The (g — 2) M constraint is 
applied up to the level of 3 a for low and 2a for high tan /? cases respectively. 


15 This is especially so for models with high tan/3 which are also in better agreement with the (g — 2)^ constraint. 
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(a) (b) 

Figure 7: Direct detection results for the LHLS (a) and LMLS (b) scenarios with tan (3 = 6 (blue) 

and 30 (cyan). Colours and conventions are the same as in Fig. 6. 

The results for the SI direct-detection are shown in Fig. 7(a) and Fig. 7(b) corresponding to the 
LHLS and LMLS scenarios for the analyses of Fig. 3 and Fig. 4 respectively. The colours (blue and 
cyan for low and high values of tan /3 as used in Figs. 3(a) and 3(b) respectively) and conventions 
used in Fig. 7(a) are the same as in Fig. 6(a). Fig. 7(a) shows only the points satisfying the major 
constraints, which in this case, lie in the range 200 GeV< m ~ 0 < 350 GeV. On the other hand, 
Fig. 7(b), which corresponds to Fig. 4(a) and Fig. 4(b), shows the above type of allowed points for 
m~o > 300 GeV. These points give rise to the LUX limit by a similar amount as before. There 
are a few points in both the figures, lying in the Higgs resonance region which are still allowed by 
all the major constraints. This region of parameter space, having very small values of cr^, will be 
fully explored by the XENON1T experiment in future. 

7.3 LHHS 

The results for the SI direct-detection are shown in Fig. 8 corresponding to a scenario of a higgsino- 
dominated xf for the LHHS analysis of Fig. 5 with tan/3 = 30. In the LHHS analysis we consider 
masses of the left and the right sleptons to be larger than m~± by 200 GeV while a choice of 
M 2 = 2/i is made to have higgsino-domination in XiVX 2 - The LSP is however bino-dominated be¬ 
cause of the choice < ^u. Fig. 8 shows only the allowed points that satisfy the main constraints. 
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Parameter points with m~o above 80 GeV satisfy these criteria. The same is true for a small region 

Xi 

around Higgs pole annihilation zones i.e. m ~o — ■^ L . cr^l for the parameter points with m ~o above 
80 GeV exceed the LUX limit but they are still within an order of magnitude. On the other hand, 
the allowed points in the Higgs resonance region lie much below the LUX limit in a region to be 
conclusively probed by the XENON1T experiment in near future. 



(a) 

Figure 8 : Direct detection results for the LHHS scenario with tan/3 = 30 for higgsino dominated 
Xi scenario. The LUX and XENON1T limits are shown as black and red lines. The points allowed 
by the WMAP/PLANCK and collider constraints are shown. The (g — 2)^ constraint is applied up 
to the level of 3cr. 


8 Spin-dependent (SD) direct detection cross-section and indirect 
detection reach for muon flux 

A larger higgsino content of the LSP, as explored in this analysis, can potentially be interesting 
for indirect detection of DM. This is principally due to a larger spin-dependent x — p scattering 
cross-section which results from a large Z\X coupling Czxx = 1-^13 ~~ -^ 14 1, where N tJ refers 
to the elements of the neutralino diagonalizing matrix. The coupling Czxx is a measure of higgsino 
asymmetry. Thus, on one hand, is already large because fi and M\ are not too far away from 
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each other. On the other hand, a^ is also in the larger side due to increased amounts of higgsinos 
within the LSP. The above enhancement in the scattering cross-sections of both the types in turn 
causes a loss of energy of a DM particle so that their velocities may go below the escape velocity. 
This results into gravitational capture of DM within the dense region of an astrophysical object. 
The captured LSPs then undergo pair annihilations. Because of an increased higgsino content 
within the LSP, the pair annihilations may lead to highly energetic neutrinos within the Sun. The 
resulting muon flux out of the charged current interactions from the neutrinos coming from the Sun 
may effectively be probed at the IceCube [39,40] experiment. In SUSY, neutralino annihilations at 
tree level would not produce neutrinos. However, the latter may arise from gauge bosons, heavy 
quarks, r-leptons etc. The neutrinos may thus have a broad energy distribution with the energy 
being limited to an appreciable fraction of the mass of the LSP. Neutrinos from bb or t + t~ are 
the primary channels when the LSP is lighter than M\y. However, due to a large threshold, the 
above neutrinos may not be suitable for detection. For massive neutralinos, LSP pair annihilation 
would produce gauge bosons, top quarks or Higgs bosons. A neutralino having a significant amount 
of higgsinos may pair annihilate to produce gauge bosons. This, in turn, may become a suitable 
source for high energy neutrinos. 

In order to estimate the capture cross-section in relation to the DM annihilation cross-section 
for the Sun, one considers the time evolution of N DM particles, 


dN 

dt 


c-c a n 2 . 


(4) 


Here C measures the rates at which DM particles are captured. C A relates to the strength of 
depletion due to DM annihilation. The annihilation rate Ta in turn is related to Ca via Ta = 
\C a N 2 [27,67,68]. 

Solution of Eq.4 results into Ta = \C A N 2 = 4Ctanh 2 (t/r) where r = 1 /yJCC A . Within the 
MSSM and for objects like the Sun (that would correspond to large annihilation and large capture 
rates) one finds that for the Solar age of t = t® = 4.5 x 10 9 years, it is justified to assume t/r » 1. 
This leads to Ta = \C, an equilibrium scenario out of capture and annihilation [69]. However, this 
is hardly possible for a less massive object like the Earth for which a captured LSP would have a 
much smaller escape velocity and where one has dominance of the spin-independent interactions 
in the DM-nuclear scattering. This leads to a weaker indirect detection signal in general [27]. In 
contrast to the above, the Sun is a massive object with a much larger escape velocity for the LSP. 


32 



At the same time, both SI and SD cross-sections are important for the capture of DM particles 
within the Sun [70,71]. Suitable models are used to relate the capture cross-section to SI and SD 
type of DM-nuclear cross-sections. Thus, a measurement of muon flux effectively sets limits on 
both SI and SD cross-sections [27, 71] . Refs. [72-74] may be seen for further details of setting the 
above limits and the associated degree of model dependence. 

8.1 LHLRS and LMLRS 

Figs. 9(a) and 9(b) show the results of scanning the aforesaid pMSSM parameter space in scatter 
plots of (Tp.® vs m ~o in the LHLRS and LMLRS models corresponding to Fig. 1 and Fig. 2 re¬ 
spectively. Blue and cyan points all of which satisfy the relic density, the collider limits and the 
(,g — 2)^ data (up to the level of 3 a for low and 2 a for high tan /3) correspond to tan/3 = 6 and 30 
respectively in each of the figures. Limits derived from the present as well as future reach of the 
IceCube experiment [39,40] are shown in black and red lines respectively. The higgsino models are 
more sensitive to the current IceCube data than the mixed models , as can be seen from the figures. 
Clearly, the final IceCube reach will exhaust the parameter space while most of the scatter points in 
general lie at most within an order of magnitude below the presently derived bound from the same 
experiment. Some points for the low tan /3 scenario in the higgsino model are tantalizingly close to 
the present limit. However, the surviving points from the Higgs resonance region for Fig. 9(b) lie 
way below the reach of even the final IceCube measurement. 

Fig. 10(a) and Fig. 10(b) show the results of muon flux 4^ in relation to m~o in the LHLRS 
and LMLRS models corresponding to Fig. 1 and Fig. 2 respectively. Blue and cyan points, all of 
which satisfy the major constraints, correspond to tan/3 = 6 and 30 respectively in each of the 
figures. Limits on the muon flux from the present IceCube data and its future reach are shown in 
black and red lines respectively. 

The scatter points in general lie within an order of magnitude below the presently derived 
bound from the IceCube. There are some points in the low tan /3 case for the higgsino model which 
lie very close to the current data. Clearly, the final IceCube reach will fully explore the parameter 
space, except the Higgs resonance region in Fig. 10(b), which is beyond the reach of even the future 
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(a) (b) 

Figure 9: Plot of spin-dependent direct detection cross-section vs the mass of the LSP for two 
cases namely (a) LHLRS : when xt is principally a charged higgsino or (b) LMLRS : when xf is an 
appreciable mixture of a charged higgsino and wino. Only the points that satisfy the relic density 
limits, the collider bounds and the (g — 2)^ data (up to a maximum of 3cr level for low and 2cr for 
high tan/3) are shown. Blue and cyan points correspond to tan/3 = 6 and 30 respectively in each 




(a) (b) 

Figure 10: Plot of muon-flux vs m~o for two cases : the (a) LHLRS and (b) LMLRS scenarios. 

Only the points that satisfy the relic density limits, the collider bounds and the (g — 2) M data (up 
to 3cr level for low and 2cr level for high tan /3) are shown. Blue and cyan points correspond to 
tan/3 = 6 and 30 respectively in each of the figures. The present and future IceCube limits are 
shown in black and red lines respectively. 
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IceCube measurement. There is a high degree of correlation between and for a neutralino 
with a significant amount of higgsino mixing. 

8.2 LHLS and LMLS 

We compute the and <3?^ in Fig. 11 and Fig. 12 for tan/3 = 6 and 30 in the LHLS and LMLS 
scenarios corresponding to the analyses of Fig. 3 and Fig. 4 respectively. The results show that the 
valid parameter space may be effectively probed in the future IceCube measurements. However, 
even the future IceCube reach can not probe the Higgs resonance regions in these models since the 
corresponding and values are too small. 



(a) (b) 

Figure 11: Plot of spin-dependent direct detection cross-section vs the LSP mass for the two cases 

: the LHLS (a) and LMLS (b) scenarios. Colours and conventions are the same as in Fig. 9. 


8.3 LHHS 

Considering the LHHS scenario associated with Fig. 5 we compute and in Fig. 13(a) and 
Fig. 13(b) with tan/3 = 30 for a higgsino dominated xf ■ Only relic density and collider satisfied 
points along with the (g — 2)^ constraint applied up to the level of 3er are shown. The present and 
future IceCube limits are shown in black and red lines respectively. Future IceCube experiment 
can fully probe the allowed parameter points. However, the h-resonance region remains beyond the 
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(a) (b) 

Figure 12: Plot of muon-flux <3?^ vs the mass of the LSP for two cases namely the LHLS (a) and 

LMLS (b) scenarios. Colours and conventions are the same as in Fig. 10. 


reach of even the future IceCube bounds. 



(a) (b) 

Figure 13: Plot of SD direct detection cross-section and muon-flux vs the mass of the LSP for 

the LHHS scenario with tan /3 = 30 for a higgsino dominated X\ ■ Only the points satisfying the 
relic density limits, collider constraints and the (g — 2)^ constraint applied up to the level of 3er are 
shown. The present and future IceCube limits are shown in black and red lines respectively. 
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9 New constraints on nig and future search prospects 

In this section we consider a scenario where the gluinos, in addition to the EW sparticles, are 
relatively light. We choose several BPs from the LHLRS (Figs. 1(a) and 1(b)), LHLS (Fig. 3(b)) 
and LHHS (Fig. 5) models. Each BP captures the essential features of the corresponding model. 
We then compute the new m-g limits in each case using the ATLAS data for dedicated squark-gluino 
searches. The published data [4-8] are based on a series of conference reports [75-78]. There are 
a few differences between the published results and the earlier analyses so far as the details are 
concerned. However, there is no qualitative difference among the two sets of exclusion contours. Our 
analysis is based on Refs. [75-78]. This exercise also enables us to anticipate the probable signatures 
of different scenarios at future colliders. The characteristics of each BP are briefly described in the 
following lines. The underlying sparticle spectra and the relevant BRs are presented in Tables. 5-7. 
We like to mention that for BP2, BP4 and BP 6 we consider tilted scenarios with the slepton masses 
driven closer to m~±. As already discussed in Secs. 4.1 and 4.2, constraints from slepton searches 
are weakened for tilted scenarios. This allows the above-mentioned BPs which represent the Higgs 
resonance regions in the corresponding figures to become acceptable parameter points. 

The point BP1 (BP2) is from the upper (lower) branch of the two regions satisfying the 
WMAP/PLANCK data in the LHLRS model (see Fig. 1(a)). From the BRs for BP1 in Table 6 
and 7, it follows that in spite of the suppressed couplings of higgsino-dominated xf, X .2 an d X 3 
with appropriate fermion -sfermion pairs of the first two generations, only the trilepton signal is 
potentially viable for the LHC Run-II. Moreover, the BR of the invisible mode X 2 —>• vv is quite 
large in this case. This is a generic feature of the branch under consideration. Thus, e + e _ —>• l/dr 
events are enhanced by XiXii X 2 X 2 an d vv final states [79]. Such events can be searched for at the 
future e + e~ colliders using the single photon tag from initial state radiation [80]. In contrast, the 
invisible decays of x% are rather suppressed in all other higgsino models studied in this paper. 

BP2 is consistent with the slepton search data with a tilted slepton mass (see the discussions 
on the LWLRS~± and LHLRS~± models in Secs. 4.1 and 4.2). Here both X 2 an d X 3 decay into Zx\ 
and hxi while xt decay into W±Xi modes with large BRs. Thus, the Whxi signal will strongly 
compete with the reduced trilepton channel during the next phase of the LHC experiments. It may 
be noted that the current LHC constraints on the former signal in the wino models [16] are rather 
weak, even if the underlying decays are assumed to occur with 100% BR. It is also worth recalling 
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BPl 

BP2 

BP3 

BP4 

BP5 

BP6 

BP7 

Taken from 

la 

la 

(Tilted) 

lb 

lb 

(Tilted) 

3b 

3b 

(Tilted) 

— 

Mi 

256 

59.5 

202 

62 

227 

62 

60 

m 2 

573 

1000 

522 

1000 

541 

960 

541 

h 

286 

500 

261 

500 

270 

480 

271 

m ?? 

229 

57 

187 

61 

208 

60 

58 

m xt 

281 

501 

258 

503 

269 

483 

268 

m x° 

295 

503 

269 

503 

282 

483 

270 

m x% 

303 

510 

272 

511 

282 

490 

282 

M- d 

(e,n) L 

266 

394 

232 

395 

249 

380 

473 

M~° 

(e,M)a 

265 

394 

232 

395 

2000 

2000 

473 

m~ 

n 

260 

387 

199 

360 

249 

380 

457 

MP 

V 1S 

255 

386 

219 

388 

236 

372 

466 


0.09 

0.11 

0.1 

0.13 

0.09 

0.12 

0.09 

a S i(pb) xlO -9 

15 

0.13 

5.4 

0.03 

7.5 

0.03 

0.2 

a s u USY xlO” 10 

A 4 

5.2 

1.8 

34 

9.6 

29 

9.2 

17 


Table 5: The sparticle spectra corresponding to different BPs chosen from Fig. 1 and 3. All masses 
are given in GeV. Only the closest figures are referred with the tilted scenarios, where the figures 
themselves represent regular scenarios with slepton masses at the midway between the masses of 
the LSP and xf- 

that this signal, in fact, does not look interesting in all the wino models except the LWHS model. 

The point BP3 (BP4), consistent with all the major constraints, is taken from the upper (lower) 
branch of the region allowed by the WMAP/PLANCK data of Fig. 1(b). The scenario represented 
by BP3 has a light t\. As a result, xt ~ X 2 pairs dominantly decay into final states involving 
multiple rs. Thus, searches with improved r-tagging may reveal this signature [59]. However, a 
small but non-negligible trilepton signal can not be ruled out a priori. Thus, final states involving 
all of the three generations of leptons in different proportions, especially at an e + e” collider, may 
be a hallmark of this model. 
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Decay Modes 

BP1 

BP2 

BP3 

BP4 

BP5 

BP6 

BP7 

9 -+ Xitt 

6.6 

7 

3.5 

6 

4 

5.4 

2.6 

-> x°R 

17 

18 

12 

14 

14 

14 

14 

->• Xlti 

13 

16 

14 

12 

12 

13 

14 

-> X°4,ti 

1.6 

- 

1.7 

- 

1.5 

- 

1.5 

-> X°ibb 

1.2 

3 

2.3 

3 

3 

2.4 

1.4 

X> 

- 

0.4 

4.4 

6.5 

5 

6.5 

5.4 

x> 

0.5 

0.3 

5 

6.3 

4 

6.2 

5 

x> 

2.7 

0.2 

2.6 

- 

2.5 

- 

2.4 

->Xitb 

44 

52 

42 

50 

44 

50 

44 

x£tb 

9.6 

- 

9.4 

- 

9.0 

- 

8 

Xi 

32 

3 

37 

18 

69 

21.4 

- 

-*■ n v T 

8 

1.5 

46 

15 

1.4 

0.4 

- 

T2V t 

- 

- 

- 

2.6 

- 

- 

- 

— >■ I'll 

53 

4 

15 

2.6 

28 

3.6 

- 

Ilvi 

6.8 

1.2 

1.8 

0.6 

2.6 

0.8 

- 


- 

90 

- 

61 

- 

74 

100 


Table 6: Dominant decay modes and BRs of gluino and chargino for different BPs taken from 
Fig. 1 and Fig. 3. 

BP4 similarly represents the generic features of the Higgs resonance region of Fig. 1(b). As in 
the previous case, final states with varieties of leptons belonging to all of the three generations is a 
feature of this scenario. However, the signal Whxi , if kinematically allowed, also looks promising 
for the LHC Run-II. Since the probability of r-rich final states is still sizable, observing both rh\i 
and Whxi could be a snroking-gun signal. The former signature, though challenging at a hadron 
collider, is expected to be rather straightforward at an e + e^ machine. 

BP5 (BP6) captures the generic features of the upper (lower) branch of the region allowed by 
the WMAP/PLANCK data of Fig. 3(b) quite well. Both the points are consistent with all the major 
constraints provided the slepton mass is tilted towards m ~±. BP5 also represents a scenario where 
all three generations of leptons may appear in the final state in different combinations, although 
a clear r-dominance is noticeable since the electroweakinos are higgsino-dominated. However, for 
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Decay Modes 

BP1 

BP2 

BP3 

BP4 

BP5 

BP6 

BP7 


2 

4 

17.6 

2.64 

0.6 

3.6 

- 

vv 

57.6 

1.2 

0.5 

0.8 

18 

1.2 

- 


10 
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22 

0.8 


- 

- 

— > t±t t 
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3 

58 

28 

82 

20 

- 
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22 
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1 

6 

- 


- 


- 

65 

- 

38 

- 

48 

62 
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- 

22 

- 

22 

- 

28 

38 

X° 3 ^lI T 

24 

- 


- 

48 


- 
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0.6 

8.4 

0.96 

6 

1.2 

- 


40 

0.4 

1.6 

0.4 


- 

- 

—> 

19.2 

- 

86 

22 

46 

20 

- 


12.8 

1.4 

4 

10 

- 

- 

- 

X°ih 


22 

- 

22 

- 

26 

20 

X°iZ 

- 

74 

- 

43 

- 

53 

80 


Table 7: Dominant decay modes and BRs of \2 an d X 3 for different BPs taken from Fig. 1 and 
Fig. 3. 

BP 6 the T-dominance is relatively mild. The Wh\i signal, as and when sufficient luminosity 
accumulates, appears to be an attractive option. 

BP7 is from Fig. 5 and is consistent with all the major constraints. Here both the trilepton and 
the Wh\i events are expected to show up. However, it follows from Table 6 and 7 that the relative 
rate of the later class of events are expected to be larger. As in the LWHS model, the LHHS model 
promises signatures other than the conventional trileptons from chargino-neutralino production. 

The gluino decay BRs for different scenarios are also presented in Table 6 . Since xt, X 2 an d X 3 
are higgsino-dominated, the final states with third generation of quarks overwhelm those containing 
leptons or light quark jets. Thus, it is expected that the search channels with tagged b-jets would 
yield the most stringent mass bounds on m~ from current data, ft also follows that these are the 
best channels for gluino search in the higgsino model during the LHC Run-11. 

Before we compute the revised gluino mass limits in the higgsino models , we summarise the 
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ATLAS SUSY searches sensitive to gluino pair production. After analysing the LHC Run-I (£ ~ 
20 fb -1 ) data, the ATLAS collaboration interpreted the results in the n-leptons + m-jets (with or 
without b tagging) + fir channel with different integral values of n and m for various simplified 
models [75-78]. For the inclusive jets + 0 1 + fir channel, depending on jet multiplicities, they 
defined five inclusive analyses channels (labelled as A to E) [75]. The selection criteria used for 
11 signal regions (SRs) are summarised in Table 1 of Ref. [75]. In the absence of any significant 
excess, an upper limit on the number of events ( Nbsm ) from any Beyond the Standard Model 
(BSM) scenario was presented for various signal regions. The most effective signal regions for our 
analysis and the corresponding upper limits on Nbsm for C = 20.3 fb" 1 are presented in Table 8. 


Channel 

Most effective 

Observed upper limits on 


signal Regions 

Nbsm (95 % CL) 


SRC-Medium 

81.2 

Jets +0/ + fir [75] 

SRD 

15.5 

SRE-Medium 

28.6 


SRE-Tight 

8.3 

Jets +11 + fir [76] 

Inclusive 6-jet(e) 

4.6 

Inclusive 6-jet(//) 

3.0 

Jets+2 SSI -31 +fir [77] 

SR3b 

3.9 

Jets (3b) +0-1 1 + fir [78] 

SR-ll-6j-B 

3.0 


Table 8: The most effective signal regions for our analysis and the corresponding upper limits 
on Nbsm at 95 % CL with C = 20.3 fb -1 in the jets + 01 + fir channel [75], jets + 11 + fir 
channel [76], jets+2 SSI — 31 + fir channel [77] and 01 +jets (3b) + fir [78] channel. 

We adopt the analysis of “hard single-lepton” (n = 1) from Ref. [76]. In this channel, the 
ATLAS collaboration defined six inclusive and six binned signal regions treating electrons and 
muons independently. Details of the signal regions are summarised in Table 4 of Ref. [76] . For our 
case, the most effective signal regions are inclusive 6-jet (electron) and 6-jet (muon) (for the upper 
limits on Nbsm in these two channels, see Table 8). 

For the same sign (SS) dilepton analysis, the ATLAS collaboration considered either two iso- 
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lated leptons (e or /i) with the same electric charge, or at least three isolated leptons (31) [77]. The 
SB,3b signal region yields the best limits for the higgsino models. In this signal region, SS or 3Z 
events are chosen with at least five jets and at least three b-jets. Corresponding upper limit on 
Nbsm for L = 20.3 fb -1 is 3.9 [77]. Details of the selection criteria for the other signal regions are 
presented in Table 1 of [77]. 

Next, we will briefly discuss the most important channel - jets (at least 3 b-jets) + 0-1Z (l = 
e, g) + Ifrr [78] which gives the most stringent bounds on m-g in the higgsino models. Selection 
criteria for the 9 signal regions are listed in Tables 1 and 2 of Ref. [78] and upper limits on Nbsm 
at 95 % confidence level (CL) are presented in Table 5 of Ref. [78] . The most effective signal region 
is SB — 11 — 6 j — C 16 , characterized by large Ifir and at least six jets which includes at least three 
b-tagged jets. 

For electron, muon and jet identification, lepton-lepton isolation, lepton-jet isolation etc., we 
follow the ATLAS prescription as described in Refs. [75-78]. For b-tagging, we use the Pt dependent 
b-tagging efficiencies presented by ATLAS collaboration in Ref. [81] . After reconstruction of objects, 
we adopt all the signal regions defined by different selection criteria, introduced in Refs. [75-78]. 
For validation purpose, we also match the number of events and efficiencies of different cuts used 
for different signal regions in Refs. [75-78] with the ATLAS results. 

Using PYTHIA (v6.428) [36] we generate the signal events in various channels from gluino pair 
production for the chosen BPs. For the NLO gg pair production cross-section calculation we use 
PROSPINO 2.1 [42] with CTEQ6.6M PDF [82]. By comparing the simulated number of events 
with the corresponding upper limits on Nbsm in the appropriate signal region, we calculate the 
new limits on rri~ for different scenarios represented by BP1 - BP7. 

The revised limits on m-g in different higgsino models are presented in Table 9. As expected, 
the strongest limits come from search channels involving tagged b-jets. The limits are practically 
independent of the choice of the slepton masses. It may be recalled that in the wino model the best 
limits come from the jets+lZ+^r channel. This table also illustrates the importance of multichannel 
search for the higgsino model. The size of jets+OZ + dftr signal can potentially distinguish some of 
the higgsino models from the others. However, our results are based on the generic strategies for 
squark-gluino searches devised by the ATLAS collaboration using tagged b-jets. In the LHC Run- 
16 signal regions are classified as A/B/C depending on fr and m e //. 
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Points 

Limit on m-g (GeV) 

jets +0/ + fir [75] 

jets +11 + pr [76] 

jets +21 + pp [77] 

0Z +jets (3b) + pr [78] 

BP1 

675 

1125 

1250 

1340 

BP2 

1080 

1175 

1135 

1360 

BP3 

815 

1100 

1180 

1320 

BP4 

1050 

1160 

1135 

1345 

BP5 

770 

1105 

1210 

1330 

BP6 

1075 

1160 

1135 

1345 

BP7 

980 

1130 

1135 

1325 


Table 9: Limits on m-g for different BPs using the ATLAS jets + 0/ + pr data [75], jets + 11 + 
pr data [76], jets +21 + pp (SSD) data [77] and 0/ +jets(3b) + p}p [78] data. 

II, more dedicated searches, e.g., the detection of a Higgs boson in a gluino decay cascade may 
provide more definite information on the underlying higgsino model. 

10 Conclusion 

The focus of this paper is on the phenomenology of the higgsino and the mixed models of the 
electroweakinos and to compare and contrast them with that of the corresponding wino models 
studied in Ref. [24]. In this concluding section we summarize our main results in the light of the 
three major constraints (see Table 10). 

To give the readers some feelings for the numerical values of the revised LHC mass limits we 
note that in the LHLRS (for low tan/3) and LHLS (for both values of tan/3) models the lower 
bounds on m~± are 380 GeV and 360 GeV respectively. These bounds are significantly weaker than 

Vl 

the similar bounds in the corresponding wino models which have the ballpark values of ~ 600 GeV. 
It is interesting to note that the entire exclusion contour in the LHLRS model with high tan /3 is 
superseded by the theoretical constraints. However, for small LSP masses, the bounds from the 
slepton searches translate into stronger bounds: m~± > 650 GeV (LHLRS) and 600 GeV (LHLS). 
As discussed in detail in the text these bounds get weaker in the tilted LHLRS~± and LHLS~± 

*1 Xi 
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Model 

Trilepton Constraint 

Processes Leading to Correct 

Relic density 

The (g — 2) M Constraint 

Satisfied at 

LHLRS 

Degraded w.r.t. the LWLRS case 

$ix ( {^w+w-,tlzz, Zh 

Z/h resonance annihilation 

ri coannihilation (only for high tan (3) 

(a) 3a (low tan /3) 

(b) < 2a (high tan/3) 

LHLS 

Degraded w.r.t. the LWLS case 

xlxi^w+w-^lzz 

Z/h resonance annihilation 

Same as above 

LHHS 

Degraded w.r.t. the LWHS case. 

xlxi^w+w-.zz.ti 

Z/h resonance annihilation 

< 2a for high tan /3 

LMLRS 

Degraded w.r.t. the LWLRS case 

but stronger than the LHLRS model. 

XiXi -+W + W-,U 

Z/h resonance annihilation 

ri coannihilation (only for high tan/3) 

(a) 3a (low tan /3) 

(b) < 2a (high tan/3) 

LMLS 

Degraded w.r.t. the LWLS case. 

but stronger than the LHLS model 

x?x? ->W + W~,tt,ZZ 

7j/\y resonance annihilation 

Same as above 


Table 10: Summary of the impact of the three major constraints on the models analyzed in this 
work. Here, LWLS corresponds to the Light Wino and Light Left Sleptons model. Similarly, 
LWLRS refers to the Light Wino and Light Left and Right Sleptons and LWHS corresponds to the 
Light Wino and Heavy Sleptons scenarios respectively. These wino models were discussed in detail 
in Ref. [24], 

models yielding m~± > 450 GeV (LHLRS) and 410 (LHLS) GeV. For higher LSP masses both the 
bounds obtained directly from the trilepton searches and those deduced from the constraints in 
the slepton sector become relaxed and eventually disappear for certain LSP masses which for each 
model can easily be read off from the figures concerned. In the LHHS model the bound is m~± > 
175 GeV for negligible LSP masses. This is rather weak even in comparison with the corresponding 
limit in the LWHS model which is the most relaxed limit among the wino models. We have also 
considered the LHC constraints in the mixed models. As expected, the LHC limits lie in between 
the corresponding ones in the wino and the higgsino models. However, no qualitatively new feature 
emerges from this analysis. 

In all models considered in this paper and in Ref [24] with low tan (3 the predictions for (g — 2)^ 
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are consistent with the data only at the level of 3 a after imposing the LHC constraints. The 
constraint can be more effective only in high tan /3 scenarios. 

In the parameter spaces of the higgsino and the mixed models there are two distinct branches 
allowed by the WMAP/PLANCK constraints, as in the wino models. In the upper branches the LSP 
pair annihilations into various channels turn out to be the dominant DM relic density producing 
mechanism as illustrated in Table 10. Only in the LHLRS and the LMLRS models with large tan (5, 
the LSP-stau coannihilation is important. In contrast, the upper branches in the wino models are 
dominated by different coannihilation processes. 

In the lower branches of the higgsino and the mixed models analyzed in this work, the DM 
production is mainly due to LSP pair annihilation via the h-resonance. The APSs are larger in 
models with tilted slepton masses for reasons explained in the text. In contrast, this mechanism is 
generically under pressure in the wino models with heavy sleptons either due to the above tension 
with the (g — 2)^ constraint or due to the LHC constraints or both. 

The inclusion of the WMAP/PLANCK and the (g — 2)^ constraints in our analysis severely 
restricts the APSs by imposing both upper and lower mass bounds in most of the models studied 
here. At high tan /3, most of the wino, mixed and the higgsino models have narrow APSs surviving 
all the major constraints. We now summarize our main findings regarding the prospects of having 
novel signatures at the LHC Run-II and at the ILC. 

In the higgsino models x 2 or X 3 decaying into Zh\\ with large BR are rather common (see 
the examples in Table 7). They occur even if the sleptons are lighter than the electroweakinos. 
In contrast, these decays occur with large BRs in the wino models with heavy sleptons provided 
the DM constraints are relaxed. The discovery of light sleptons together with the observation 
of the Whfrp events due to chargino-neutralino production during the LHC Run-II could be the 
hallmark of the higgsino models. If we focus on this signal in parameter spaces consistent with the 
(g — 2)^ and the DM relic density constraints, then only zones with high tan /3 and DM relic density 
production via LSP pair annihilation into the h-resonance are acceptable. On the other hand this 
DM producing mechanism is generically disfavoured in the wino models with high tan /3 as noted 
earlier. 

In some regions of the APSs, especially in the upper branches of the regions allowed by the 
WMAP/PLANCK data in the m~± — m~ 0 plane, the conventional trilepton channel appears to be 

Xi Xi 
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the best bet in many models considered in this paper and in Ref. [24]. With improved r-tagging 
efficiencies, the final states with multiple r’s as analysed in Ref. [59] for models with high tan /3 
may provide alternative/complementary search channels during the LHC Run-II. Several novel 
signatures which could be relevant at the ILC have also been discussed in Sec. 9. In particular, the 
comparison of the constraints obtained in this paper and the ones in Ref. [24] clearly indicates that if 
the ILC indeed operates at around 500 GeV during its first run as planned, then the electroweakinos 
in the higgsino model have larger probabilities of being within its striking range (see, for e.g., Fig. 5). 

Assuming the gluinos to be light in addition to the electroweak sparticles while all squarks 
are heavy, we have revisited the gluino mass limits in Sec. 9 using the ATLAS data. As in the 
wino model , this analysis emphasizes the importance of multichannel searches. It follows that the 
conventional jets + 01 + Ifir signal has the poorest sensitivity in a wide variety of the higgsino 
models (see the results in Table 9 for different BPs in Sec. 9), whereas jets + 11 or 21 + 1/irp signal 
has a better sensitivity. The best channel for probing the higgsino models involves multiple tagged 
b-jets. The gluino mass limits are stronger in general than the ones obtained in the corresponding 
wino models. This observation may help to formulate the future strategies for gluino searches in 
the context of the higgsino models and distinguish between the wino and the higgsino models, if a 
signal is seen. Moreover, one can distinguish among various higgsino models consistent with the 
major constraints by the relative rates of jets + 0 l + l/dr and jets (3b) + 0 l + l/dr events (see Table 
8 ). 

We next summarize the prospects of direct and indirect detection of DM in the context of the 
higgsino models. There is a significant bino-higgsino mixing in the LSP in the higgsino models 
especially in the portions of the parameter spaces consistent with the DM relic density constraint 
and characterized by relatively small m~± — m~ o. The spin-independent direct detection cross- 
section a^ p is larger in all such cases compared to the corresponding wino models. Although the 
cross-sections exceed the LUX limits in most of the cases, they stay within an order of magnitude of 
the same limits. However, the cross-section for the points from the Higgs resonance region allowed 
by the major constraints (e.g. in Figs. 6(b), 7 and 8) lie well below the LUX limit. All the models 
can be probed by the future XENON1T experiment irrespective of the degree of uncertainties in 

u? 1 . 

XV 

The spin-dependent cross-section a?9 and muon flux values for the neutrino signals in the 
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higgsino models are enhanced compared to their wino model counterparts. Most of the points in 
almost all the scenarios are allowed by the present IceCube data. In some cases the situation seems 
very interesting since the values of <7-^ and muon flux lie very close to the present experimental 
bound (see e.g. Figs. 9(a), 10(a), 11(a), 12(a) and 13). All the other cases would decisively be 
probed by the future IceCube searches. For the points representing the /i-resonance region in Figs. 
9(b), 10(b), 11, 12 and 13 the values of these observables are too small to be detected even by the 
future IceCube reach. 
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